Kinematics of the swimming of Spiroplasma 



On 
O 
O 



Jing Yang,^ Charles W. Wolgemuth/' ^ and Greg Huber^' ^' ^'B 

^ Center for Cell Analysis & Modeling, University of Connecticut Health Center, Farmington, Connecticut 06030, USA 
^Department of Cell Biology, University of Connecticut Health Center, Farmington, Connecticut 06030, USA 
^Department of Mathematics, University of Connecticut, Storrs, Connecticut 06269, USA 

(Dated: June 29, 2009) 

Spiroplasma swimming is studied with a simple model based on resistive-force theory. Specifically, 
we consider a bacterium shaped in the form of a helix that propagates traveling- wave distortions 
which flip the handedness of the helical cell body. We treat cell length, pitch angle, kink velocity, 
and distance between kinks as parameters and calculate the swimming velocity that arises due to 
the distortions. We find that, for a fixed pitch angle, scaling collapses the swimming velocity (and 
the swimming efficiency) to a universal curve that depends only on the ratio of the distance between 
kinks to the cell length. Simultaneously optimizing the swimming efficiency with respect to inter- 
kink length and pitch angle, we find that the optimal pitch angle is 35.5° and the optimal inter-kink 
length ratio is 0.338, values in good agreement with experimental observations. 
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Helices are profoundly involved in the swimming of 
many types of bacteria. For example, the classical pic- 
ture of bacterial swimming involves long helical filaments 
extending out from the inner membrane of the cell 
driving the bacterium forward by their rotation. But 
helices are also used in other more ingenious ways. In 
Treponema pallidum^ the spirochetal bacterium responsi- 
ble for syphilis, rotation of helical filaments encased be- 
tween the bacterial cell wall and the outer membrane 
leads to periodic undulations that drive this cell through 
water and other viscous fluids. (This is not, presum- 
ably, how it crossed the Atlantic Ocean [2].) The Lep- 
tospiraceae rely on a helical body plan and rotation of 
tightly-coiled fllaments to drive their motility [3]. Gen- 
erally speaking, it is the dynamics and shape together, 
the "swimming strategy" , that is the key to understand- 
ing the motility of helical bacteria. For the flagella-less 
family of plant pathogens represented by Spiroplasma^ 
the precise propulsive strategy was shrouded in mystery 
until careful observations [4, 5] revealed that these tiny 
organisms processively flip the chirality of regions of their 
helical bodies, generating pairs of kinks moving down the 
cell body (Fig. [H^a)). The internal propulsive mecha- 
nism underlying this peculiar swimming strategy is still 
unclear, as is the quantiflcation of its motility. 

Spiroplasma swimming is realized by the propagation 
of a pair of kinks along the body axis of the cell, and 
the bacterium is propelled by the hydrodynamic force as 
the fluid associated with the kinks moves rearward with 
the "body wave" . Kinks start at the same end of the cell 
(front), and travel toward the other end (back). As the 
kink propagates, the cell changes direction through an 
angle related to the pitch angle of the helix. Because of 
the unbalanced viscous drag between different portions 
of the cell body separated by kinks, eventually the cell 
swims in a zigzag path. The helical transformation re- 
quires a cell to rotate about its body axis during kink 
propagation [1, 0, 0] • 

In this Letter we use the experimentally-observed ge- 




FIG. 1: (a) Spiroplasma in differential interference contrast 
(Image courtesy of J. W. Shaevitz). (b) Schematic of helical 
geometry. D is the helix diameter, D -\- d is the outer coil 
diameter (D = 0.138/im, d = 0.20/im). Here the bend angle 
a and pitch angle ip satisfy a = n — 2ip. 



ometry of Spiroplasma to explore its swimming behavior 
for a range of kinematic deformations. We consider a 
helically-shaped cell of flxed pitch angle and helix diam- 
eter. To swim, the bacterium flips the handedness of 
its helix beginning at one end. This change in helicity 
propagates toward the back of the cell with a uniform 
velocity, Uk = 10.5 ± 0.3 /im/s [5]. After a time t^, the 
front of the cell reverts to its original handedness and this 
change in helicity also propagates toward the back of the 
cell with the same kink velocity. The distance between 
the kinks is Lk = Uktk, and only an integer number of 
turns is considered. We treat Uk^ Lk, and the cell length, 
Lc, as parameters and use resistive- force theory (RFT) to 
calculate the swimming speed and trajectory as a func- 
tion of these parameters. Our model is complementary 
to the work by Wada and Netz Q where slender-body 
hydrodynamics [9] was used to explore the dependence 
of swimming speed on the pitch angle, tp. By varying the 
distance between kinks, we examine both the geometry 
and kinematics of the swimming process. It was found 
experimentally that the kink velocity Uk 0, IH and the 
inter-kink distance I//c [1] show a remarkable consistency 
- what determines this preference? Our results demon- 
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strate that the maximum propulsive efficiency is achieved 
when the pitch angle is around 35.5° (in agreement with 
0]) and the ratio between inter-kink length and cell-body 
length is around 0.338, which is consistent with most ex- 
perimental observations and suggests that the kinematics 
of Spiroplasma has been evolutionarily tuned. 

Any helix can be described by its pitch, p, and he- 
lix radius, r. One helical repeat in the arclength s 
is i = \/ pP" ^ Air'^r'^. When two helices of differing 
pitch are concatenated, such that the tangent and nor- 
mal vectors are continuous, a kink is formed [6, 
The bend angle a between the two axes in this con- 
struction is a = TT — 2iIj (Fig. [I{b)). The parame- 
ters employed in the calculations are p = 0.62 /im, 
r = D/2 = 0.069 /im, and = 35° 0, For the 
numerical integration, the cell body is discretized into 
straight-line elements. According to RFT [III, [l2[, the 
resistance force acting on a unit length of rod is given by 
f = - [27r/i/ln(2/e)] (21 - AA)-U + [ln-^(2/e)] , where 
U is the local velocity, I is the unit tensor, e is the slender- 
ness ratio of rod diameter to length, A is the unit tangent 
vector, and /i is viscosity. The total force and momentum 
on the unrestrained cell body is zero, i.e., J^"" i{s)ds = 0, 

f^"" 'M.{s)ds = 0. The translational and rotational veloc- 
ities of the cell body can be solved for from the resulting 
system of linear equations. To verify that the implemen- 
tation of the model has been properly made, the trans- 
lational velocity of a rigid helix with angular velocity co 
around its body axis is calculated with the same numer- 
ical code. Since the total force on the rotating helix is 
zero, i.e., i{s)ds = (L is the total length of the rigid 
helix), the translational velocity of the helix is U = Bico^ 
where B = sin^ cos?/;/ [27r(l + sin^ t/^)] The nu- 

merical results match this resistive-force solution with 
relative errors less than 0.05% under our typical spatial 
discretization. A slender-body theory (SBT) calculation, 
using the straight line elements, was also performed to 
account for hydrodynamic interactions. Based on the 
aforementioned model, a number of dynamical questions 
on the swimming of Spiroplasma can be formulated and 
addressed. 

Typical experimental and simulated displacement 
curves during one swimming cycle of Spiroplasma are 
shown in Fig. [2](a). Here the displacement is defined 
as the projected center-of-mass position. For compar- 
ison with the experimental data, our numerical simula- 
tion results are also smoothed with a 175-ms boxcar filter 
Using the parameters given above, our simula- 
tions agree well with the corresponding experimental re- 
sults [l4], giving an average velocity of about 3 /im/s, 
close to the observed cell velocity, 3.3±0.2 /im/s 

It is observed experimentally that the speed is de- 
creased when two kinks are present as compared with 
when only one kink is present. This key feature of the 
displacement curve is reproduced in the numerical sim- 
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FIG. 2: (a) A simulated displacement curve (Lfc=2^, Lc=7^, 
and [/fc=12-£/s). (b) Velocity ratio and the corresponding en- 
velope angle, as a function of the length ratio for different 
body lengths, from RFT (fit by solid line); Maximum veloc- 
ity ratio is located where envelope and pitch angle are equal. 



ulations (Fig. [2l^a)). Our simulations show that the cell 
velocity increases linearly with kink velocity, and that 
the ratio of cell velocity to kink velocity is constant for 
fixed inter-kink length, with the proportionality constant 
varying with inter-kink length, consistent with the exper- 
imental reports 0, IH- 

In order to delineate the effect of dynamically changing 
geometry on Spiroplasma motility, two dimensionless pa- 
rameters are introduced. One is the ratio of cell velocity 
Uc to kink velocity Uk] the other is the ratio of inter- 
kink length Lk to cell body length Lc. We find that all 
the data fall on a single curve after the results from dif- 
ferent cell-body lengths and inter-kink lengths are taken 
together (Fig. [2](b)). The curve describes a velocity ra- 
tio that increases with length ratio at small values, and 
then begins to decrease with length ratio after reaching 
a peak value. From Fig. [2](b), the maximum velocity 
ratio of the data-fitted curve is around 0.323 when the 
length ratio is about 0.32. The fact that the velocity 
ratios and length ratios observed in the experiment fall 
close to the aforementioned value suggests that wildtype 
Spiroplasmas have optimized their geometry to achieve 
the "fastest state" in the course of evolution. The study 
of Wada and Netz [8] suggests an optimized pitch angle 
as well, the treatment of which goes beyond the reach of 
RFT for pitch angles greater than roughly 35° [l5|, 
The full optimization over both length ratio and pitch an- 
gle (Fig.[3|) using our SBT code gives a maximum velocity 
ratio at L^/Lc = 0.338 and t/; = 35.5°. These values are 



consistent with the experimental values 0.303±0.051 an 
34.8±1.4° respectively 
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FIG. 3: Velocity ratio as a function of kink length ratio and 
pitch angle. Contour lines are velocity- ratio isolines. Experi- 
mental observations are centered within the error ellipse. 

The velocity ratio can be calculated in RFT in the 
limit that the length ratio approaches one. As the length 
ratio goes to one, there is effectively only one kink prop- 
agating along the cell body at any given time. In this 
case, Spiroplasma can be treated as a filament consisting 
of two helical sections, one right handed and the other 
left handed. The two helical sections rotate around their 
respective axes with the opposite sense. The total trans- 
lational velocity is the vector sum of these two compo- 
nents with an interior angle 2tp. From geometrical con- 
siderations, the kink velocity Uk can be represented as a 
function of rotational velocity cj, i.e., Uk = 2u;r/ tdiUip. 
Using the previously-mentioned translation velocity of a 
rigid helix {B£u) the total velocity can be estimated to 
be 2B£ujcosiIj, and the corresponding velocity ratio re- 
duces to 271 B = sin^ cos (l + sin^?/;). For a pitch 
angle of 35°, this ratio is around 0.202, not far from the 
simulation result of 0.24. At a length ratio of one, the 
kinematic asymmetry between the left-handed inter-kink 
length and the right-handed extra-kink length vanishes, 
and a chirality variable x = (^c — Lk)/{Lc -\- Lk) can be 
used to study the velocity ratio on both sides of x = 
[12] . We note in passing that the entire velocity-ratio 
curve is in excellent empirical agreement with a quartic 
polynomial pj| 



Uc/Uk = 0.24 + 0.73(/)^ - l.eic^"^ 



(1) 



(the solid curve in Fig. [2] (b)) where (j) = tan~^ x- 

Not only translation, but also rotation of the cell body 
must satisfy the dynamic balance, i.e., zero force and 
torque. At the beginning of a kink cycle, one end of the 
cell-body axis traces out a circular arc (see Fig. [4|). The 
maximum angle swept out, we term the envelope angle. 
We use this concept in Fig.[2](b), where the dashed curve 
shows that the envelope angle increases with length ratio. 
Fig. EJb) also shows that the velocity ratio achieves its 




FIG. 4: The trace swept out by the cell-body axis at eight 
different inter- kink lengths 2-9 i (for cell length 10 i). Three 
snapshots of the cell body are superposed in the leftmost trace 
(Lk = 2i). Yellow lines indicate the trajectories of the two 
kinks. The envelope angle is indicated. 

Experiments on Spiroplasma do not observe more than 
two kinks along the cell body at one time. It is there- 
fore interesting to consider the swimming efficiency as 
a function of the number of kinks per cycle in order to 
determine whether a two-kink cycle is optimal. A stan- 
dard definition of swimming efficiency compares the en- 
ergy dissipated by dragging the appropriate rigid cylin- 
der through fluid to the power dissipated by the actual 
swimmer [8|, |12|, ll9|. This definition compares a swim- 
ming object to a non-swimmer and also does not account 
for the energy dissipated to create the swimming motion. 
Since an unknown mechanism is responsible for produc- 
ing the traveling kinks in Spiroplasma^ we assume that 
there is a fixed net energy, SE, required to produce a pair 
of kinks [lO] and that this mechanism sets a uniform kink 
velocity. Under these assumptions, we hypothesize that 
biology tries to maximize the distance traveled per unit 
energy expended. Therefore, we consider the fuel mileage 
7 = distance traveled per energy consumed {nSE + the 
energy dissipated by the fluid, where n is the number 
of kink pairs per cycle), which we compute in terms of 
nanometers traveled per ATP burned First, ignor- 
ing the energy to create the kinks {SE 0), we calculate 
7 as a function of the length ratio (solid curve in Fig. [5]). 
The fuel mileage from different cell lengths can all be 
collapsed onto one curve by constructing a scaling fac- 
tor from a power of the slenderness ratio e (= d/Lc) as 
e-^j{Lk,Lc) with ^ :^ 2. In both RFT and SBT, we 
find that the fuel mileage is maximized at a length ratio 
of around 0.3, although the methods different velocities 
lead to shifted 7 values: 7rft/7sbt — 3 — 3.4. Using 
RFT, we find that a single double kink cycle requires 
50-100 ATPs worth of energy and the cell moves about 
26 nm/ATP, whereas SBT gives 150-340 ATPs per cycle 
and the cell moves 8 nm/ATP. Interestingly, the latter is 
the fuel mileage of a single kinesin molecule [22|. When 
6E is small, we find that the fuel mileage for a 4-kink cy- 
cle is larger than for a 2-kink cycle. However, when 5E 
is greater than 40 ATPs, then a 2-kink cycle is optimal. 
Therefore, we expect that the mechanism producing the 



4 



double kink must require at least 40 ATPs per cycle. 
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FIG. 5: Fuel mileage as a function of the length ratio. The 
solid curve and circle data points show the case with Lc = lOi 
(left ordinate). Using the scaling factor e~^, the results from 
all cell lengths can be collapsed. Here e is the slenderness 
ratio and the exponent 6 = 1.90 ±0.01 [23| . 

In this Letter, we investigated via a numerical study 
the swimming of Spiroplasma with a model based on 
resistive-force theory and the kinematic deformations of 
Spiroplasma. The key features of Spiroplasma swimming 
are well reproduced with this model, and the simulation 
results agree quantitatively with the observed experimen- 
tal data. The important effects of dynamic geometry on 
the motility of Spiroplamsa are explored, and the motility 
efficiency is shown to assume the highest possible value 
given the kinematics observed in experiments. One impli- 
cation is the optimization of the geometry and kinematics 
of Spiroplasma during the course of evolution. 
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